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Abstract 

The excitation of wake waves by a relativistic homogeneous elec- 
tron bunch passing through cold magnetized plasma at equilibrium is 
studied for arbitrary values of the ratio of bunch density to plasma 
density. The analysis is based on the assumption that the magnetic 
field is sufficiently strong {ub ^ ujp, where cUp and ub are the elec- 
tron plasma and cyclotron frequences respectively). The periodic and 
nonperiodic solutions for the momentum of plasma electrons inside 
and outside the bunch are analyzed. It was shown that the presence 
of strong external magnetic field may increase the amplitude of wake 
waves and the maximum transformer ratio, the latter being reached 
at densities of bunch lower than that in the absence of magnetic field. 
The optimum conditions for obtaining the maximum values of the 
wake field amplitude and of the transformer ratio are found. 
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1 INTRODUCTION 



The plasma wake- field accelerator (PWFA), proposed by Chen, Dawson, Huff, 
and Katsouleas^ in 1985, and discussed in a slightly different context by Fainberg^ 
as early as 1956, uses electrostatic fields in a plasma wave driven by a relativistic 
electron beam. Charged particles can then be accelerated in this electrostatic 
wake to ultrahigh energies. The first experimental observation, at the Argonne 
National Laboratory Advanced Accelerator Test Facility, of acceleration by the 
PWFA mechanism has been reported in Ref.^ The one-dimensional theory of the 
nonlinear regime of the PWFA developed by Rosenzweig^"^, Ruth et al.^, Am- 
atuni et al.^'^ and Khachatryan^" predict certain advantages over the linear regime 
(for more details see reviews^^'^^). It has been shown that the transformer ratio, 
the ratio of the maximum decelerating field inside the driving electron bunch 
to the maximum accelerating field found in the wake of the driver, is enhanced 
by driving the plasma waves with an ultrarelativistic electron bunch of density 
one-half that of the plasma electrons. This interesting method of obtaining large 
transformer ratios is asserted to be more straightforward than alternatives in the 
linear regime, and also to lessen the multiple scattering that the accelerating par- 
ticles undergo in a PWFA, as the plasma used need not be as dense as the linear 
regime requires. The simplicity of the nonlinear scheme for obtaining high trans- 
former ratios, as well as the inherently high accelerating gradients, also makes 
this method interesting as a candidate for high-energy cosmic-ray-acceleration 
mechanism. 

The high transformer ratios obtained in this nonlinear scheme depend on 
driving the plasma electron density waves to extremely large amplitudes. In 
contrast to the linear waves, the maximum energy the electrons reach during 
the oscillation is many times their rest mass energy. The electron oscillation 
period is an increasing function of this maximum energy. As the amplitude of 
the wave grows, the wave steepens dramatically, and the positive excursion in 
density becomes very large in amplitude and narrow in time. 

An important result of the nonlinear theory is the proof that the wave breaking 
limit is £^max = {mui,ujp/e) [2(75 — 1)]^^^ and it is reached at rib/rio < 1/(2 + 7^^^), 
where Up = (47rnoe^/m)^^^ is the plasma frequency of electrons, Ub and no are the 
densities of the bunch and plasma electrons, 76 = (1 — P^) is the relativistic 
factor of the bunch, = Ub/c, uj, is the velocity of the bunch. The Dawson^^ 
wave breaking limit is equal to E'max ~ 2muujp/e when 7^ w 1 (/3b <^ 1). 

The one-dimensional relativistic strong waves can be excited in the plasma 
by wide relativistic bunches of charged particles or intense laser pulses^^'^^ (when 
kpttQ ^ 1, where kp — ujp/ubi Oq are the characteristic transverse sizes of bunches 
or pulses). 

The excitation of nonlinear one-dimensional wake waves in underdense and 
overdense plasma by a relativistic electron bunch was studied in Ref.^ 
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In the present work an analogous problem has been treated in the presence of 
strong external magnetic field Bq. A qualitatively novel effect of dependence of 
the amplitude of wake wave excited by an one- dimensional bunch on the angle a 
of magnetic field orientation with respect to the direction of bunch motion was 
found. When a is increased from zero to A = arcsin(l/76), then wake waves 
with larger amplitudes as well as larger transformer ratios are obtained. Here, 
the maximum value of the transformer ratio is attained at lower densities of the 
bunch in comparison with the case of Bq = 0. 



2 BASIC EQUATIONS 

Our treatment will start with equations of cold relativistic nonlinear fiuid theory. 
We assume that the plasma ions form an immobile neutralizing background of 
density no and define the plasma electron density n, velocity v = /3c, and the 

1 /2 

plasma frequency ujp = (47rnoe^/m) . We now write the fiuid equations for 
plasma electrons in the presence of an ultrarelativistic beam with homogeneous 
density n^, length d, moving with the velocity U;, = /3^c along the z axis and 
external constant magnetic field Bq {Bq^ — Bosina, B^y = 0, Bq^ — B^cosa, 
where a is the angle between Bq and Ub): 

1 (9E 

V X B = - — - AnenfS-Ane^f^rib, (1) 
1 flYi 

VxE=— — , V-B = 0, (2) 

COT 

V • E = -Ane {n - Uq) - Anerib, (3) 

^ + c (/3V) p = -— [E + /3 X (B + Bo)] , (4) 

OT mc 

dn , , 

_ + V • (nv) = 0, (5) 

where p = (3/ {1 — 13"^) e is the absolute value of electron charge. 

If we assume that the wave motion is a function only of the variable ^ — z — ut 
(the wave ansatz, with taken as the wave phase velocity) we obtain nonlinear 
differential equations for the dimensionless momentum and field components: 



^(0=^0—^== , (6) 
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d^ 

dpxiO 
dC 

dpyiO 



mccOp 1 



Pz 



+ 



rib 



e Xp ypf^^l + p^ - no 



e ^ cos a 



e 



Xb p,y^TT7-P^' 

1 Pz sin a — Px cos a 
f3by^l + 7^' 



(8) 

(9) 
(10) 



p - p^ 



dpz 
di 



P^ + lT {pxEx + PyEy) 



sm a 



-Py, 



(11) 

where 7^ ^ = 1 — (3^, Xp = c/up, Xb = c/ujb, ojb = eBo/mc. The components of 
induced magnetic field are determined by means of relations: E^ — /3bBy, Ey = 
PbBx: Bz — 0. Note that as it follows from Eqs. (6)-(ll), px = Py = E^ = Ey ~ 
in the absence of external magnetic field (Bq = 0) and Eqs. (8) and (11) pass 
into the well known equations derived in Refs.[9] and [15]. So, in an anisotropic 
magnetized plasma, besides the longitudinal motion, the one- dimensional bunch 
excites the motion of plasma electrons in the direction normal to that of bunch 
travel. When a — (the magnetic field is directed along u^) the magnetic field 
has no infiuence on the excitation of wake waves and we arrive at the known 
equations of motion that describe only the longitudinal motion of electrons. 

To investigate the set of Eqs. (6)-(ll) it is convenient to exclude the elec- 
tric field and obtain an equation containing the components of dimensionless 
momentum: 



de 



cos a\/l + — Pb {Px sin a + pz cos a) 



^ 1 PhbPx sill a- Pz cos a _ cos a rib 

Ap PbVTT7-Pz ~ K ^o' 

Then the component Py is eliminated from Eqs. (9) and (11). The following 
equation is obtained: 



+ 



— (Pj; sin a cos a) (13) 
di 

g (l3b\/l + p^tga + Px- pz^ga] Ex + pyEy + I3b\/l + p'^E^ 
= — ^cosa-^^ , . 

In what follows we shall assume the external magnetic field to be compara- 
tively strong so that the condition ujb'^ ojp {or Xb ^ Xp) were observed. Hence 
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we obtain the following restriction on unperturbed concentration of plasma elec- 
trons: 

Bo>3x (14) 

where Hq is measured in cm^^, Bq in kG. The condition (14) is always observed 
in the range of parameters Uq < lO^^cm^'^ and Bq < 100 kG. In the zeroth-order 
approximation in one obtains from Eqs. (7) and (9) Py — Ey — B^ — 0, 
from Eq. (10) p^sino; = p^cosa (or — pcosa, p^ = psina), i.e., the plasma 
electrons move only in the direction of external magnetic field. From Eqs. (12) 
and (13) we obtain in case of As —>■ 0: 



— (cos ay/1 + p2 - p^pj + (15) 



7^ sin^ a — I p cos a rib 

/56v^l + p2_ pcosa A2 no' 

/3i,y/l + p^ — p cos a dp 



yTTp^ di mc 



{ExSm.a-\- Ezcosa) . (16) 



To express the components of electric field E^ and E^ through the dimen- 
sionless momentum of electrons p, one is to obtain the second equation that 
connects the components of the field. One can derive it from Eqs. (4) and (5) by 
eliminating the momentum of electrons. As a result one has 

d , ^9 9^ \ mcUr, ^9 o sin auh , 
- {E, cos a + (5W,E, sin a) = '-(^lll-. (17) 

Integrating the system of Eqs. (16) and (17) taking into account the boundary 
conditions at the bunch fronts ^ — d {E^id) — Ez{d) — 0) and ^ = {Ey.{—Qi) — 
-E'a;(+0), Ez{—Q) = Ez{+0)) one obtains the following expressions for the field 
components E^ and E^ inside (0 < ^ < d, ri;, 7^ 0) and outside (^ < 0, rih — 0) 
the bunch: 

0<^<d; 



EM) - -^r, f$(e)cosa + «sin^a^^ ) , (18) 



e 1 — 7b sin a\ uq X. 



p 



m) = = ^^Jtt^ ( t(e) + cca'^i^) . (19) 

e 1 - 7ft sm^ a \ no Xp J 



e<0; 



E,{0 = $ C cosa + Z^ft^^in^a^- , 20 

e 1 — 7^sm q;\ rioXpJ 
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e 1 — 7^ sin^ a \ no X, 



where 



m = \ 



(22) 



Note that from the continuity of components of the field and momentum of 
electrons on the boundaries of the bunch there follows the continuity of function 
$(^) at { = and ^ = rf. 

Thus, Eqs.(15) and (18)- (22) determine the wake wave in the magnetized 
plasma for large value of the external magnetic field. For a = we obtain from 
Eqs. (18)-(22) that = By = 0, and Eqs. (15), (18) and (20) coincide with 
equations given in Refs. [9, 15], that describe longitudinal oscillations of the 
plasma electrons. When a has values different from zero, the nature of solutions 
of Eqs. (15) and (18)-(22) is drastically changed. First, as we shall see below, the 
behavior of solution of Eq. (15) depends on the fact whether the angle a is greater 
or less than a definite value A = arcsin(l/7b). Second, owing to an anisotropy of 
plasma due to the presence of an external magnetic field, besides the longitudinal 
components there arise transverse components (as well as an induced magnetic 
field). Furthermore, owing to the anisotropy of plasma the components of the 
wake field (Eqs. (18)-(22)) contain the Coulomb field of the bunch (the second 
terms in Eqs. (18)-(22)) that depends along with other parameters also on the 
angle a. 

Below we shall analyze the Eq. (15) inside and behind the bunch; inside, on 
and outside the anisotropy cone a — A. 



3 ELECTRON MOMENTUM INSIDE THE BUNCH 
FOR ARBITRARY VALUES OF THE RA- 
TIO UB/no 



Now let us consider the equation of plasma motion (15) inside the bunch. In this 
range < ^ < 0?, nb — const 7^ 0. The integration of Eq. (15) with due regard 
for boundary conditions p{d) — $(d) = will give 




(23) 



where 



A — 1 — sin^ a 




(24) 
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In Eq. (23) one is to take the sign plus (minus) when p(^) decreases (increases) 
with increasing As it follows from the Eq. (23), for solution of Eq. (15) it 
is necessary that the radicand in (23) assume only the positive values. Hence 
follows a restriction on the momentum of plasma electrons: p < 0. Taking into 
account this restriction one is to take the sign " minus" during the integration of 
Eq. (23). Making allowance for boundary conditions at the front boundary of 
the bunch we obtain from (22) and (23) 

[ Pb\/l + 7f — r] cos a] drj 

; ^ =, (25) 

^/l + rj^^A - A^l + rf - Brj 

that determines an implicit dependence of p on ^. Note that the expression 
obtained in [9] will follow from (25) at a = 0. 

To proceed with an analysis of the expression (25) one has to know whether 
a is greater or less than A and the dependence of a on the density of bunch. 

First consider the case when a is inside the anisotropy cone (a < A). 




3.1 Inside the Cone: a < A 

Consider first the case of small densities of the bunch. 



3.1.1 The Case of ^ < 7 ^°^""^" = h 

no '0 cos a ^ 

In this case < B < A and one can easily show that the range of variation of p 
is determined by the relation — po < P < 0, where 



2(3b£ ^cosa 
l-Pie^' l-72sin2«-f^cos2a- ^^^^ 

Then the dependence of p on ^ has a periodic behavior with oscillation amplitude 
Po- 

The integration of (25) will give an implicit dependence of p on ^ for afore- 
mentioned values of the magnetic field orientation angle a and density of bunch: 



rf-f _ 1 /2(l + ;cosa))V2,Ml + -u«,)^^^_^j_ 



Xp 1 - Y 1 - 72 sin^ a 1 VI + Pb£ 

- {(3b + cos a) \Jl- Pb^P - a/I + > 



where 



V'^arcsinf ^yW^7^T7) , ^ ^ J < 1, (28) 
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E{ip,k) is the elliptic integral of the second kind, < (3be < 1. When p — po, 
assuming in Eq. (27) that ^ = one obtains the length of bunch do or the wave 
half length inside the bunch 



Ain — 2(io 



(1 + e cos a) 



1 - /56£ V (1 - 7' sin' a) (1 + Pbs) 



E{k), 



(29) 



where E{k) is a complete elliptic integral of the second kind. 

The dependence of wavelength inside the bunch on the bunch density is shown 
in Fig.l for values of a = (sohd line), a — 0.2 A (dashed hne), a — 0.4 A (dotted 
line) and < rib/no < bi. 

For values of ratio rih/nQ > foithe only constraint on the region of allowed 
values of p is the condition p < 0. Therefore, the dependence of p on ^ will get 
nonperiodic and p{^) will increase infinitely with ^. 

3.1.2 The Case of ^ = 61 

no ^ 

In this case A — B and integration of expression (25) gives 



(cos a + Pb 



2]j 2(3b{l-'ylsm^a) 

|^i-p-v/TT7(Vi + p'-p-2 



(30) 



cos a ~ f3b 
cos a + Pb 



+ 



+ ln 



One can see from Eq. (30) that at large values of c? — ,^ ^ Ap, i.e., with the 
distance from the front boundary of the bunch |p(^)| grows linearly depending 
on d — ^. 

3.1.3 The Case of 61 < ^ < 72cosa+^ ^ ^ 

^ no cosa ^ 

In this case —B < A < B. Taking into account this inequahty on obtains from 
Eq. (25) 



A„ 



^"^^ !^{B + \A\) {cos a + Pb)\fA^^A^AT^^^p+ (31) 



52 



+ {B cos a + l3b\A\) 



B-\A\ 

~7W 



F(x,«;)-2Vi?i?(x,«:) 



where 
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X = arccos 



B+\A\ 
2B 



F{x, k) is an elliptic integral of the first kind , k < 1. 
3.1.4 The Case of ^ = 62 

no ^ 

In this case A = —B. One obtains from the general form of Eq. (25) 



— = —1= { cos a + /3b - - (cos a- Pbji -x/T+p^ + p\ 
p \J1B \ \_ 2 \ / 



X 



(32) 



(33) 



X \j -sjl^ (P- - p- \ - ^ (cos a - f3b) arctg y^yr+ p2 _ p _ i| . 



3.1.5 The Case of ^ > bo 



no 



Here A < —B. The calculation of the integral in Eq. (25) gives 



A. 



A^-B^ 



X 



(34) 



X 



\A\ cos a — PbB + {(3b\A\ — 5 cos a) 



_ 52 _ A vTTp2 - Bp 



B^Jl + p^ + Ap 



+ 



+ {B cos a + (3bA) 



\l \/A^ - 52 ^AE{a, r) + 



Vv^2 - 52 



F(a,r) 



where 



(T = arcsm \ \ , , r 



2V^2 _ 52 

V>1^- 52-^4' 



(35) 



In all cases under consideration in items 3-5, at large values of d — ^ ^ Ap 
|p(^)| is a quadratic function of the distance to the front boundary of the bunch 
d — ^. At a = the expressions obtained in items 1-5 coincide with the results 
given in Ref. [9]. 



3.2 On the Cone: a = A 

In case of a = A the Eq. (15) for the momentum of electrons is simplified ant 
takes the following form: 
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de 



no' 



(36) 



The integration of Eq. (36) with due regard for boundary conditions gives the 
following expression for p(^): 



2no \ \ J 


rib I 
2n() \ 




2 


nb_ ( ±z£' 
2no \ Ap ^ 


2 

+ 1 





(37) 



As it follows from Eq. (37) for arbitrary densities of the bunch the momentum 
of electrons increases monotonically with the distance from the front boundary 
of the bunch. 



3.3 Outside the Cone: a > A 

At large values of angle a {a > A) the value of A is always negative. As follows 
from an analysis of Eq. (23), in such a case p is a monotonic function of ^ and is 
described by Eqs. (31)-(35) for different values of the bunch density. So, when the 

conditions rih/uQ < 621 n^b/i^o = ^2 and Uh/no > 62 are approached, the implicit 
dependences of p on ^ are given by Eqs. (34), (33) and (31) respectively. 

4 THE FIELDS INSIDE THE BUNCH FOR 
ARBITRARY VALUES OF THE RATIO m/no 

In this Section the Eqs. (18) and (19) describing the induced electromagnetic 
fields inside the bunch were analyzed. Consider at first the case a 7^ A. We 
shall begin with an analysis of a transverse component of electric field E^. In the 
previous Section it was shown that the momentum of electrons inside the bunch 
meets the condition p(^) < 0. As follows from this inequality, Eq. (7) and the 
boundary condition E^^d) = 0, inside the bunch E^ is positive and monotonically 
increases with distance from the front boundary of the bunch. The function Er^ 
(as well as By) takes on the maximum value at the rear boundary of the bunch, 
at ^ = 0. Prom Eq. (18) one obtains for the maximum value of E^: 



^.max = ^-^^^^V^-^ ±2i/yyl-/10T7(O)-5pO + — ^cosa . 

e 1 — % sm a \ ^ no Ap / 

(38) 

In Eq. (38) one takes the sign "plus" ("minus") if p(^) is decreased (increased) 
with increasing of ^. In particular, if n^/no < 61 (the momentum of electrons 
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inside the bunch is periodically changing with ^ and is given by Eq. (27)) and 
d — Ai„/2, then p(0) = — po and the maximum value of as obtained from (38) 
is 

_ mcujp rib Ain p'hf sin n cos a 
e no 2Xp 1 — % sm a 

In all other cases discussed in the previous Section, E^^jmax is determined by 
the Eq. (38) with the sign "minus". 

The longitudinal component E^ of the field satisfies the Eq. (8). In case of 
nb/no > cos a /{(5b + cos a) = be, dE^/d^ < and Ez{i) is positive and mono- 
tonically increase with the distance from the front boundary of the bunch. In 
case of nb/no < K the longitudinal field E^ first monotonically increases with the 
distance from the front boundary of the bunch, than reaching the maximum at 
p(^) — —p, where 

(3,9 g= V^o ^ (40) 



monotonically decreases to the rear boundary (^ = 0). 

The maximum value of longitudinal field E^ in the point ^ follows from ex- 
pression (18), where the substitutions p — — p and ^ — ^ are made: 



= , """T"^" { cosaJA - A^lTf + Bp - (3^1 sin^ 

l-% sm^ a y ^ no Xp 

'(41) 

where the value of {d — ^)/Xp is expressed through p by means of expressions (27), 
(30), (31), (33) and (34). The value of longitudinal field at the rear boundary of 
the bunch (^ = 0) follows from the Eq. (18): 

£^.(0) = ^ -"^^pI^ f±2V2 cos aJA - A^l + p\Q) -~B^) + sin^ a"^^ 
l-% sm^ a \ ^ noXp 

(42) 

In the general case the maximum value of longitudinal field inside the bunch 

is 



Ezmax — max 



E,{0)\; E, . (43) 



In particular, if Ub/uQ < hi and d = Ain/2 the value of E^ is determined by 
the Eq. (41), and for E^(0) we obtain from (18), (23), (26) and (29): 

^^(0) = A, 
^ ' e l-j^sm^ano2Xp ^ ' 
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The case cos a = Pf, (or 7b sin a = 1) is to be considered separately, as in this 
case the Eq. (18), (19) make no sense. Indeed, it follows from expressions (22) 
and (36) that both the denominator and the determinant of the expressions (18) 
and (19) tend to zero when a ^ A. The electromagnetic field inside the bunch 
in this case may be found from Eqs. (7), (8) (or Eqs. (7) and (17)) and (37). 

As a result we obtain 



A. 



n-b ( <tzl\ ' 
2no \ \ J 



+ 1 



{rib/no) 



arctan 



nb d- ^ 
2no Xp 



(45) 



^^(0 - 7 



1 mcUr, 



d - 4 no V no ^) \ X, I ^ ^ 



4ni, 
no 



At; 



{rib/ no) 



arctan 



nb ( ±<\ 
2no \ J 



Ub d — ^ 

2no Xp 



(46) 



+ 1 



The values p, ^ and are determine by means of the following expressions 

e 



X„ 



- 1 



no 



(47) 



(48) 



1 mcujri 



4 e y (nb/no) 



\ 



I _ 2n6 
no 



lU/l-^-2(l-^))+(49) 
no no 



+ arctan 



I _ 2n6 

no 



Below we shall consider the equation of motion of plasma (15) outside the 
bunch (C < 0). 
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5 THE WAKE FIELD 



5.1 Inside the Cone: a < A 

To find the wake field Ex{C) and £"^(0 behind the bunch < 0), one has to 
integrate Eq. (15) for n?, = taking into account the condition of continuity of 
fields and momentum p at the rear boundary of the bunch = 0. In this case 
the fields are determined by means of expressions (20)- (22), where now 



Here 



A - Dp{0) 
1 — 7^ sin^ a 



+ Vi + P^(o)(i- 



.4 



1 — % sin a 



- 1 



(51) 



is the maximum allowable value of the momentum p{C,)- As is seen from expression 
(50), the momentum p(^) (hence, the fields Ex{^), -E^(0 and function $(^)) is 
a periodic function and changes in the range —pmax < p(0 ^ Pmax- The wake 
wavelength is determined from expressions (22), (50) and has the following form 



A 



out 



\/r^^o7sii?a 



X _ p2 



E{p) 



Kip) 



where 



(52) 



'1 



1 



P 



(53) 



\/l + P^ax + l' 

The dependence of wake wavelength (Aout) on its amplitude pmax is shown in 
Fig. 2. It is seen that the wavelength increases with the amplitude and angle a. 

As it follows from Eqs. (20)- (22) and (50), (51) the transverse {Ex{C)) and 
longitudinal {Ez{C)) wake waves change within the ranges Ei^ < E^i^) < S2x, 
^iz < ^ziO < ^2z, respectively, where 



^2x 



mcUp Pl'-^l sin a 
e 1 — 7^ sin^ a 



2nh 
no 



cos a 



1/2 



Ub d 
cos a — -— =F 
no Xp 



(54) 



cos a (^1+7(0) - l) - MO) 
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^2z 



moujp/e 



l — ^l sin^ a 



q2 2 ■ z 

P676 Sin a 



rib d 
no Xp 



± 



(55) 



2nb 
no 



cos^ a 



1/2 



^cosa(^A+7M-l) -/?6P(0) 



It is seen from obtained expressions that \Six\ < E^^x and > \£2z\- Thus, 
the maximum values of transverse and longitudinal fields behind the bunch are 
£2x and \£iz\ respectively. Here the value of transverse field varies with respect to 
some positive quantity (the first term in Eq. (54)), and the value of longitudinal 
field varies with respect to some negative quantity (the first term in Eq. (55)). 

Now consider the efficiency of wake field excitation for different values of 
the density and thickness of the bunch. If the momentum of plasma electrons 
inside the bunch is a periodic function (the bunch density meets the condition 
of subsection Al of Section III, Ub/no < &i), then inside the bunch the variation 
range of p is defined by the relation — po < P < 0, where po is determined 
by expression (26). In case when an integer number of wavelengths is interposed 
within the bunch {d = (n+l)Ain, where n = 0, 1, 2 ...), then it follows from results 
of subsection Al of the Section III that p(0) = 0, and in virtue of expression (51) 
Pmax = 0. In this case (as in the absence of magnetic field [9]) the momentum 
of electrons for < is zero, and the density of plasma is coincident with the 
equilibrium value riQ. As follows from expressions (20), (21), and (50) and (51), 
the components of electric fields for ^ < arc constant and are determined by 
the second terms of expressions (20) and (21). Thus, when the condition of 
d — {n + l)Ain is observed in an anisotropic plasma, a constant electric field 
may arise behind the bunch that is perpendicular to the external magnetic field 
according to expressions (54) and (55). 

In case oid = (n+l/2)Ain, the amplitude of wake field oscillations is maximum 
Pmax = |p(0)| = Po- The maximum values of electric field components are 



^(out) 



mcuj. 



2 cos a 



_y/(l-7^in^a) (1-/3,^) 



+ 



(n + l/2)Ain l3bJb sin^ a 
Xp cosa (1 + ecosa) 

(56) 



p(out) ^ ;^'^P^/32_2 
a; max 



£/?676 sina 



Ain n + 1/2 



y(l-7,^sin^a) (1-^,V)^ A,l + £COsa-^ 



(57) 

The condition of positivity of plasma density n(^) (Eq. (6)) imposes a restric- 
tion on possible values of the momentum of plasma electrons at ^ < 
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OO < p(0 < Pmax < 



Pblb 



(58) 



\/l - lb 

that, in its turn, limits the values of electron momentum at the rear boundary of 
the bunch = 0). For example, for the value of bunch density rib/riQ < hi Eq. 
(58) yields 



Pmax(O) 



< p(0) < 0. 



(59) 



For a bunch, the density of which satisfies the condition nf,/no < 6i, the 
electron momentum inside the bunch is a periodical function. In this case from 
expressions (26) and (59) we obtain a restriction on the density of bunch, that is 
more rigid than the condition rib/no < bi: 



Ub ^ 

no ~ 



1 — 7^ sin^ a 



cos a ^2 cos a + \J%'^ — sin^ (x^ 



= 6o < h- 



(60) 



For the density of bunch rib/riQ = bo, the maximum values of fields have the 
following forms: 



^(out) 



mca;„ 



7.^-1 



V2 



cosa 



y (l - 7^ sin^ a) cos a + ^/T 



Pblb sin^ " (Ain/Ap) (n + 1/2) 



cos a ^2 cos a + 



2 • 2 
— sm a 



7^ sin^ a 



mojj, 



^/^b 76 sin a 



^/2|3b 



(1-7^ sin^ a) ^76 cos a + ^/l 



482) 



7^ sin^ a 



^ {n + 1/2) 

276 cos a + a/I — 7^ sin^ a 

The dependence of wake waves on ^ is shown in Figs. 3-6 for the values 
7b = 60, d = Ain/2, Ub/no — bo{a) (Figs. 3 and 4) and Ub/no — 0. 860(0;) (Figs. 
5 and 6). Note that in Fig. 3 the curve calculated for the value a — Q coincides 
with the result obtained in Ref. [9]. 
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In case of even higher densities of the bunch rib/no > bo, there is also a hmit 
on the thickness of the bunch d < rfmaxi where d^^x is determined from expression 
(25) if in the latter we put ^ = and p — — Pmax(O): 



\ /-PmaxCo) {(3b\/l +1]'^ + 7] COS a) drj 

rfma. = ^ / ^ , ^ (63) 

V 2 Jo .Ji^^A - A^flT^ + Br] 

When 6o < ^fe/^o < ^i; the thickness of the bunch should satisfy the condition 
< (i < (iinax or nAin — rfjnax ^ d < nAin + (iinax) whcrc u = 1, 2, 3, ... Note that 
in this case dmax < Ain/2. For rib /no > &i the thickness of the bunch satisfies only 
the condition < d < cimax- 

In Fig. 7 the dependence of maximum thickness of the bunch its density is 
shown at rib /no > bo for a = (solid line), a — 0.4 A (dashed line), a = 0.8 A 
(dotted line). As was expected, dmax decreases with increasing density of the 
bunch. 

Consider now the variation range of the value of p(0) for n^/no > bo- From 
expressions (51) and (52) one has — Pmax(O) < p(0) < 0, where Pmax(O) is now a 
function of the density of beam: 



.(0) - 



1 — % sin a 
where 



cos ay^ _ (^^ 2 _ ^-^^2 _ 



(64) 



a/1 — " sin^ n ( " cos n — \/l — " sin^ n 

G = cosa-\ — - — r —. (65) 

[nb/no) cos a 

One can find the values of electric fields at the rear boundary of the bunch 
(^ = 0) for d = dmax from expressions (18), (19), and (23), if one makes a 
substitution p = — Pmax(O) in the latter, where Pmax(O) is determined by the 
expression (64). As a result we have for the value of $(0) 

$(0) = -V2 1 7b cos q; •^/l - 7^ sin^ a - 7^ ^Gcosa - /St-^ - {%'^ - sin^ a)^ | 

(66) 

The maximum values of the longitudinal and transverse electric fields inside 
the bunch are determined by expressions (41), (43) and (66). 

5.2 On the Cone: a = A 

In case of o; = A and = the Eq. (36) takes the following form: 

d^ 



^^2vv- . P^-P)=0. (67) 
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The integration of Eq. (67) with regard to the boundary conditions gives the 
following expression for the momentum of electrons: 



rib d 



d/2-i 



nt d rf/2-g _|_ 2 
no Ap Ap 



2no Ap 



rib d _|_ 2 

TlQ Ap Ap 



(68) 



The components of electric field are determined from expressions (7), (8) and 



(68) 



no Ap 



no Ap A. 



_d_ 

Ap 



2no I A ' ^ 



"6 
no 



+ 2 



A,] 



(69) 



£;,(0 = £^.(o) + 



mcUp ^ 



rib _d_ 

no Ap 



_d_ 

Xn 



rib ( d/\' 
2no \Apy 



+ 2 



e 2Ap!ikA£VHii + 1 

no Ap Ap 



nt 
no 



+ 2 



1 

At, 



(70) 



where Ex{^) and £"^(0) are the values of electric fields at the rear boundary of 
the bunch, and can be determined from expressions (45) and (46). As follows 
from obtained Eqs. (68)-(70), the modulo of momentum and electric field grow 
linearly with the distance from the rear boundary of the bunch. 

5.3 Outside the Cone: a > A 

If the orientation angle of the magnetic field is outside the anisotropy cone, then 
the dependence of the momentum of electrons and of electric fields on coordinates 
is monotone. 

Introducing the notations 



D 
C 



7^ sin a — 1 > 



(71) 



COS a 



cosa Vl + P'(0)-1 -MO) 



-D 



and integrating the equations of motion we obtain an implicit dependence of the 
momentum of electrons on ^: 
When C> D > 0; 



(72) 
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+^^^i^^l + ^[F{iy{p),T)-F{iyo,T)+E{iyo,T)-E{iy{p),T) 

c 



^___^[F(,(p),r)-F(,o,r)] + 



Vl + p2(0) + l 



\/rT 



where 



pip) 



C-D 
C + D 



, ^{p) — arcsin 



C + D 



, i/o = i^(p(0)), (73) 



arcsin 



-P 



\/rT7+i 

When -D<C <D; 



, po^p{p{0)). 



1 



where 



cos a 



D 



-2 



{f(|-.(p),.)-f(|-.o,^)- 
£;(|-Kp),^)-^(|-^.o,^)] + 



(74) 



2 1 



(75) 



One can obtain the electric fields from expressions (20) and (21), where $ is 
determined from the expression: 



(76) 



Taking into account the fact that in the space — oo < ^ < the function p(^) 
is monotonically increasing, the sign chosen in the expression (76) is "minus". 
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6 THE TRANSFORMER RATIO 



In this Section we shall consider the transformer ratio R that is defined as the 
ratio of maximum value of longitudinal wake field to that of decelerating field 
inside the bunch R — -E'i°m2x/-E'z max [4-9]. Since the periodic wake field is excited 
only in case of a < A, below we shall consider the transformer ratio precisely for 
these values of angle a. 

We would briefly remind the results obtained for R in the absence of magnetic 
field. As follows from the results of Refs. [4-9], the transformer ratio monoton- 
ically increases with density of bunch, the maximum value i?max — \/2 (76 — 1) 
being achieved at Uh/uQ ~ l/(2-|- 1/7?,). In the linear case (n^/no <^ 1) the func- 
tion R assumes the value R = 2. When rih/rio > 1/(2 + 1/7^), the transformer 
ratio is monotonically decreasing and tends to unity, i? ~ 1, at very large values of 
bunch density. So, in this case for excitation of wake waves the optimal is a bunch 
with length d < dmax = Ain/2 and density 1/(2 -|- I/76) < rib/rio < 1/(1-1- Pb): 
for which the amplitude of wake waves and the transformer ratio take maximum 
values. 

Consider now the influence of external strong magnetic fleld on the trans- 
former ratio. In case of bunches of low density, when n^/no < bo, the maxi- 
mum value of longitudinal wake field is determined from Eq. (55), where d — 
(n-|- l/2)Ain, (n = 0, 1, 2, ...) and p(0) = — po- One can determine the maximum 
decelerating fleld inside the bunch from expressions (40) and (41). For very small 
values of the density of bunch {ub/no &i) and n = we find from expressions 
(40), (41) and (55): 



, , 2 cos^ a + uPk'Jk sin^ a 

Ro{a) — 



1—7^ sin^ a 
2 cos^ a 



(77) 



-y/(l - 7^ sin^ a) [1 + (7^ - 2) sin^ a] - 2/?^7^ sin^ a arcsin 

that for two limiting values of the angle a takes on the following form: 

{2; a^O 
3(2+.)(l-,-^)^^^ ■ ■ 

The maximum wake field as well as the transformer ratio R grow monoton- 
ically with increasing density of bunch, the maximum value of the field being 
achieved at Ub/no = bo and being determined by Eq. (61). 

As was mentioned above, at Ub/rio > bo the thickness of the bunch is restricted 
by the value dmax (see Eq. (63)). If the density of bunch is in the interval 
^'o < nb/no < be, then one can obtain the maximum value of decelerating field 
from expressions (40) and (41), where the expression {d—^)/\p itself is determined 
by Eq. (25) with p = —p. For large densities of the bunch Ub/rio > be the 
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longitudinal field inside the bunch monotonically decrease from a definite value 
at the first boundary to zero at the second boundary of the bunch. Then the 
maximum value of decelerating field inside the bunch is equal to the value of 
field at the first boundary (when ^ = 0) and is determined by expression (42). 
In the range of rib/no > bo as a whole, the maximum value of longitudinal wake 
field is \£iz\ (see Eq. (55)), and, hence, in the last two cases (Eq. (42) and Eq. 
(55)) instead of po and d one is to put their maximum allowable values (see Eqs. 
(63)-(65)). 

In the limit of large densities of the bunch (n^/no ^ be) the following expres- 
sion is obtained from Eq. (63): 



^"^^ " KA^ \l ~ " (t^cosc - 0^^^^;^) , (79) 

that enables one to find asymptotic form of the transformer ratio for large den- 
sities of the bunch: 

l-% sm^ a 

As follows from expressions (77), (78) and (80), in the presence of external 
strong magnetic field in plasma, the transformer ratio increases and for a < A 
much exceeds the value of at = 0. 

In Figs. 8-10 the dependence of R on the density of bunch [ub/no) (Figs. 8, 9) 
and on the angle a (Fig. 10) is shown for % = 60. In Fig. 8 the dotted, dashed 
and solid lines correspond to the following values of the angle a respectively: 
a — 0, a — 0.2A and a — 0.4A. In Fig. 9 the dotted and solid fines have 
been plotted for a = 0.6A and a = 0.8A respectively. In Fig. 10 the solid line 
corresponds to Ub/no = 0.7, the dashed and dotted lines respectively to Ub/rio = 1 
and to Ub/rio = 1.3. 

One can assure that the dotted line in Fig. 8, for which a — 0, coincides with 
an analogous result obtained in Ref. [9] in the absence of external magnetic field. 

Thus, it follows from the above analysis and from Figs. 8-10 that the presence 
of strong external magnetic field increases the amplitude of wake waves excited 
by the one-dimensional bunch as well as the transformer ratio. This effect is 
most strongly manifested in the case when the angle of magnetic field orientation 
is close to A (a < A). Note also that the transformer ratio is maximum at 
/no ^ bo, and with increasing a the position of this value is shifted to lower 
densities of the bunch. Hence, in the presence of strong magnetic field the larger 
values of transformer ratio may be obtained at much lower densities of the bunch 
than in the case of Bq — 0. 
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7 CONCLUSIONS 



In conclusion we shall give a brief summary of the results and discuss possible 
restrictions on the vahdity of conclusions. 

The behavior of one-dimensional electromagnetic waves excited in plasma in 
the presence of external strong magnetic field, strongly depends on the orientation 
of bunch motion with respect to the direction of magnetic field. At large values of 
angle a (a > A = arcsin(l/7b)) the density and momentum of plasma electrons 
, as well as the excited electromagnetic fields vary in the space monotonically. 
Here the electromagnetic field increases linearly with the distance from the bunch. 
Obviously, this result is valid only in the one-dimensional model and in case of an 
allowance for the transverse size of the bunch (i.e., in the three-dimensional case) 
the fields will tend to zero at large distances. In Ref. [16] an excitation of three- 
dimensional linear wake waves was considered in the magnetized plasma. It was 
shown here that at large distances from the driving bunch the fields nevertheless 
tend to zero although slowly (by the power- law dependence). So, the linear 
increase of one-dimensional wake fields behind the bunch is limited by the distance 
that is proportional to the cross-size of the bunch. 

The allowance for the transverse size of the bunch results also in the following 
additional restrictions [17]: 7^A^ <S c? and <S a^, where a is the transverse 
size of the bunch and r is the transverse coordinate. 

At small values of the angle o; (« < A) a periodic wake fields is excited 
behind the bunch, the amplitude of which increases when a ^ A. In this case 
the transformer ratio also increases, the position of its maximum being displaced 
to lower densities of the bunch. 

It follows from the aforesaid that optimal conditions for the excitation of wake 
waves are d < Ain/2 and n^tjuQ ~ &o(Qi), where h^^OL) is determined by expression 
(60). 

Another feature due to the presence of the magnetic field is the excitation 
of transverse electric and magnetic fields, that are absent in the one-dimensional 
case at So = [4-10]. 

The allowance for thermal motion of plasma electrons leads to inessential 
corrections to wake wave characteristics at temperatures practically achievable 
under laboratory conditions [18]. 
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FIGURE CAPTIONS 



FIGURE 1: The dependence of the Ai^/Ap on n^/no (0 < ni^/riQ < hi) for 
7b = 50, q; = (sohd Hne), a = 0.2 A (dashed hne) and a = 0.4A (dotted Hne). 

FIGURE 2: The dependence of the Aout/Ap on pmax for 7^ = 50, a = (sohd 
line), a = 0.6 A (dashed line) and a = 0.8A (dotted line). 

FIGURE 3: The dependence of the function eEz{^) /mcUp on^/Xp for nij/riQ = 
60(a), d = Ain/2, 7b = 60, q; = (solid line), a = 0.2A (dashed line) and a = 0.3A 
(dotted line). 

FIGURE 4: The dependence of the function eEx{C)/f^cu}p on^/Xp for Ub/no — 
bo{a), d = Ain/2, 7^ = 60, a = O.lA, a = 0.15A and a = 0.2A. 

FIGURE 5: The dependence of the function eEz{C,)/iTLCujp on$^/Xp for riiy/riQ = 
0.860(a), d = Ain/2, 7b = 60, a = (solid line), a = 0.2A (dashed line) and 
a = 0.4 A (dotted line). 

FIGURE 6: The dependence of the function eEx{^) / maop on^/Xp for Ub/nQ — 
0.860(0;), d = Ain/2, 7b = 60, q; = O.lA (solid line), a — 0.15A (dashed line) and 
a = 0.2 A (dotted line). 

FIGURE 7: The dependence of the (imax/ Ap on rih/riQ (rib/no > 60) for lb = 60, 
a = (solid line), a = 0.4 A (dashed line) and a = 0.8 A (dotted line). 

FIGURE 8: The dependence of the transformer ratio R on rib/no for 7b = 60, 
a — (sohd line), a — 0.2 A (dashed line) and a — 0.4 A (dotted hne). 

FIGURE 9: The dependence of the transformer ratio R on rib/no for 7b = 60, 
a = 0.6 A (dotted line) and a = 0.8 A (solid line). 

FIGURE 10: The dependence of the transformer ratio R on anisotropy angle 
a for 7b = 60, rib/no = 0.7 (solid line), rib/no = 1 (dashed line) and rib/no = 1.3 
(dotted line). 
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